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Abstract 

We consider the stochastic Stokes' drift of a flexible dumbbell. The dumbbell 
consists of two isotropic Brownian particles connected by a linear spring with zero 
natural length, and is advected by a sinusoidal wave. We find an asymptotic approx- 
imation for the Stokes' drift in the limit of a weak wave, and find good agreement 
with the results of a Monte Carlo simulation. We show that it is possible to use 
this effect to sort particles by their flexibility even when all the particles have the 
same diffusivity. 

1 Introduction 

Stochastic Stokes' drift is the modification of classical Stokes' drift [1] by random forcing, 
for example Brownian motion[21 E]. This process can be used for particle sorting^ E] as, 
to leading order in the strength of the wave motion, the drift from waves with different 
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spatial or temporal frequencies are additive. This is because the Stokes' drift first appears 
at second order in the wave strength and the cross-terms average to zero. Thus it is 
possible to make particles of different diffusivities drift in different directions. There are 
also oceanographic applications pE]. Some exact results are known in one-dimension p2, Ej. 

Processes of this kind, however, are not restricted to the direct fluid mechanical appli- 
cations, and occur in many other systems, for example in molecular motors, see [7] and 
the references therein. 

In this study, we consider a simple model of a flexible particle, namely an elastic 
dumbbell, for which analytical results are still possible. Such a dumbbell model might 
be used as a crude model of a flexible polymer in solution. We show that particles of 
different flexibilities (i.e. different spring constants) have different Stokes' drift velocities 
even if they all have the same diffusivity. 

Furthermore, there exist parameters for which the dependence of the Stokes' drift on 
the flexibility is non-monotonic, even for particles having the same diffusivities. So, even 
in one-dimension, particle selection is possible. 

Because, to leading order in the Stokes' drift, waves of different spatial or temporal 
frequencies do not interact, the different drifts for particles having different flexibilities 
but the same diffusivities can be made to 'fanout' (i.e. drift in different directions). This 
extends the possibilities for using stochastic Stokes' drift in particle sorting. 



2 The stochastic dumbbell process 

Consider the motion of two identical isotropic Brownian point particles X and Y, with no 
hydrodynamic interactions between them, and joined by a linear spring, with relaxation 
time A -1 , and of zero natural length. We henceforth refer to A as the spring constant 



(even though it is not in its usual units). 

The absence of hydrodynamic interactions would not be expected to qualitatively 
change the results, and including them would lead to a very much harder mathematical 
problem. Also, if hydrodynamic interactions were included, one would probably want to 
consider a more realistic model of a flexible polymer than the dumbbell model considered 
here. This, however, is left for future work. 

In the absence of a wave, i.e. in a fluid at rest, we set the diffusivity of the centre of 
the two particles to D = \<y 2 . Thus, the stochastic dumbbell equations, in Ito form, are 

dX t = (ef(X t , t) - \\{X t - Y t )) dt + 2 1 ' 2 <jdB u 

(1) 

dY t = (ef(Y u t) + \\{X t - Y t )) dt + 2 1 / 2 adW t , 
where B and W are independent standard Brownian motions, / represents the wave 
motion, and e is a dimensionless parameter, which will be used as a small parameter in 
the asymptotic analysis to follow. 

This system can be studied in several different ways. We adopt a rational asymptotic 
approach, as often used in fluid mechanics, albeit using stochastic calculus, rather than, 
say, centre manifold theory. The structure of the calculation is simple, as it uses well- 
known results throughout, and is very similar to the standard argument for classical 
Stokes' drift. 

Suppose that the time t £ [— T, T], where T = e~^T , and T is the largest natural 
time scale in the system. These natural time scales are the spring relaxation time A -1 , 
the period of the wave forcing, and the typical time for a particle to diffuse a wavelength. 

Now consider two new processes given by U t = ^(X t — Y t ), and V t = ^(X t + Y t ), where 
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'=' means denned to be equal. From (CQ), we find 

dU t = \e{f{X u t) - f(Xt, t))dt - XU t dt + adB t , 
dV t = \e(f{X u t) + f(Y t , t))dt + adW t , 
where B and W are independent standard Brownian motions. 
We now consider a formal asymptotic expansion in e for 

X t = Xi 0) + eX[ l) + e 2 xf ) + • • ■ , Y t = Y t {0) + eY t m + e 2 Y^ 2) + 

U t = Ul 0) + eU® + e 2 uf ] + • • • , V t = + eV^ + eV t (2) + 
in the limit e — > 0. We now take each order in e one at a time. 



(2) 



(3) 



2.1 Order £° 

At this order, becomes 

duf ] = -XU^dt + adB t , dV t {0) = adW t . (4) 

This implies 

where B is a standard Brownian motion, with B(0) = 0, independent of W. Note that, for 
simplicity, we have started the stable Ornstein-Uhlenbeck process with its stationary 
law; though this choice does not affect the long-time behaviour of interest, it does avoid 
starting transients. So to leading order, the centre of the dumbbell moves as a Brownian 
motion with diffusivity D = |a 2 , and explains the strange choice constants in (OQ). 
Dropping the bars, we find the following solution of ([T]) to order e°: 

X t (0) = aW t + e- xt B{\\- W At ), 

(6) 

Y t {0) = aW t - e-^B^X-W^). 
Note, however, that B and W are not the same Brownian motions as in equation (CQ). 
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2.2 Order e l 

From ([2]), at order e 1 , we find 

T(1) =-\uM + l(f(xW,t)-f(Y t W,t)), 

(7) 



(// -\uM + l(f(xl \t)-f(Y}°\t)), 

dt l(f(X?\t) + f(Yl°\t)), 

from which the corresponding expressions for X^ 1 ' and immediately follow. We, 
however, need only X' 1 ', which is given by 



(i) 



eft 



where 



(9) 



e -A Q t _ a) _ f( y t % t - «)) da, 

where taking the upper limit of the integral to be oo is equivalent to ignoring starting 
transients. Even though strictly the asymptotic series is valid only on a finite interval 
[— T, T] , extending the integral beyond T introduces only exponentially small terms that 
do not contribute to the asymptotic series at any finite order in e. Also, we are interested 
only in the Stokes' drift to leading order, and this appears at order e 2 , as it does in the 
classical Stokes' drift problem. 

2.3 Order e 2 

From (121), at order e 2 , we find 



dt 



d ^-\{ax?\t)xv+aY?\t)YV } , 



(10) 



dt 

where /' is the derivative of / with respect to its first argument. From (flOl) . the corre- 
sponding expressions for X^ and Y^ immediately follow. 
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3 The stochastic Stokes' drift of a dumbbell 



We now consider a special case in which 



f(x, t) = u cos(fcc — ujt + 



(11) 



where u, k and u are constants, and is a random phase, uniformly distributed on [0, 2tt). 
The Stokes' drift will not depend on 0, but this choice of phase avoids starting transients. 
In this section, we are interested only in finding Stokes' drift 



V = lim t~ L (V t - Vq) . 

t—+oo 



(12) 



From (HI) and (T71), we see that there is no contribution to the Stokes' drift at orders s° 



and e 1 , so we find 



V = e 2 E 



(2) 



dt 



(13) 



where we have used the Ergodic Theorem to replace the time limit in f|T2l by an expec- 
tation, and this is why it was important to avoid starting transients. Using ( TTOl) . reduces 



EI) to 



V = e 2 E 



l(nX?\t)xM+nY{°\t)YM) 



+ 



(14) 



Note that both of the terms in the expectation of ( 1141) have the same law, thus 



V = e 2 E 



(15) 



From (jSJ) and ©, and writing V = e 2 V^ + ■ ■ • , we find 

E [f'(Xi°\t) [f(Xi%, t-fi)- XUl%)\ df3. (16) 

Again, extending the integral to oo removes starting transients, and the contribution 
from (T, oo) is exponentially small. To evaluate this, we need 



E 



f(xl°\t)f(xl%,t-p) 



-WkE 



sin (k{x[ 0) - Xf\) - uo(3 



(17) 



as the other terms in the complex exponentials average to zero. As x[^ — X t _« is Gaussian 
and has zero mean, this reduces to 



E 



f(xP,t)f(xl%,t-(3)] = ± M 2 fcsin(^)exp (-§A; 2 £[(Xf - X^%) 2 }) , (18 



where, from 



E[(XP - X[%f\ = a 3 (^ + A- 1 (l-e-^)) 



(19) 



To complete the evaluation of ffTol) . we use © to show 



exp(— Act) 



(20) 



x £ 



da. 



The first term in the expectation on the right-hand side is as in (1171) but with different 
times, so we are essentially left with finding 



K( a + (3) = E /'(Xf, t)f(Y t Z^, t-a-(3) 



(21) 



which itself is only slightly different from (|17[) . As X^ — Y^_ T is Gaussian with zero 
mean, we find 



K(t) = lu 2 ksm(cor) exp (-^[(Xf - Y®)*\ , , 



(22) 



where r = a + (3. From 



E[(x (o ) _ y W )2] = a 2 ( r + A -i(l + e-^)) 



(23) 



We are now ready to evaluate 



M((3) ee E [f'(X?\t) (f(xl%,t- 0) - Xu£\ 



M{(3) 
\u 2 k 



sin(o;/3) exp (-|JfeV ((3 + A" 1 (l - e~ A/3 ))) 

A exp (-Ac* - \k 2 a 2 (a + [3 + A" 1 )) sin(cu(a + (3)) sinh ^ 



(24) 



from (fTol) . Gathering together what we have just evaluated in (|18p and (|22p . we find 



k 2 a 2 
2A 



(25) 



Thus ( fT5l) and ( TTBT) imply that the Stokes' drift to leading order in e is 

poo 

V = e 2 M(/3)dp + 0(e A ). (26) 
Jo 

Note that in the case of sinusoidal forcing the first correction is 0{e i ) rather than 0(e 3 ) 
as replacing e by — e is the same as a phase shift, so does not change the Stokes' drift. In 
fact, all terms with odd powers of e are zero, in this case. 



4 Strong and weak spring limits 

We briefly consider the large and small A limits. In the large A limit, the particles coalesce, 
and we expect to recover the point particle result[2]. As \/u — > oo, with k 2 a 2 /u fixed, 



M -» \u 2 k sin(a;/3) exp(-±&V/3), (27) 



which reduces (EBP to 



On the other hand, small A means that the spring is weak and the particles can 
be separated by many wavelengths. If the dumbbell is being used to model a flexible 
polymer, this is not a very physical limit, but we mention it here as this dumbbell model 
may have applications outside of fluid mechanics[2]. 

In the limit \/uj — > 0, with k 2 a 2 /u fixed, we find 

M -> lu 2 ksm(top)exp(-k 2 a 2 P), (29) 

which reduces ( 1261) to 

i 9 u 2 kui , . 

V = ^FTO + -- (30) 

Note that this has the same form as the single particle result but with different constants. 
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5 Numerical results 

We already know that we can separate particles with different diffusivities using stochastic 
Stokes' drift [2J, and we could always choose wave parameters to make the flexibility of 
the particle unimportant. So we present here numerical results for particles that all have 
the same diffusivity to highlight the scope for sorting particles where only the flexibilities 
are different. 

A graph of the leading-order asymptotic result for the Stokes' drift as a function of 
the spring strength (A), with all other parameters set to 1, is compared to a Monte Carlo 
simulation using a stochastic Euler integration with time steps of 0.001 and taken up to 
t = 10 s (see figure [T]). The agreement is good even when the 'small' expansion parameter 
e = 0.5. Other parameter values were also tested in the same way with equally good 
agreement. 

Note that in sorting applications, the characteristics of the wave motion would be, in 
general, under user control. 

6 Conclusions 

We have shown that it is possible to have a Stokes' drift that is strongly dependent on 
the spring constant of the dumbbell particles even when all the particle diffusivities are 
the same. Thus, in principle, it is possible to use this effect to sort particles solely by 
their flexibilities. 

The leading-order asymptotic result for the Stokes' drift of an elastic dumbbell agrees 
well with simulations even for quite large values of the 'small' expansion parameter. 
Furthermore, the dependence of the Stokes' drift is not necessarily monotonic in the 



9 



0.095 




Figure 1: Stokes' drift against A for e = 0.5 and all other parameters equal to 1. The 
curve is the leading-order asymptotic result and the data points are from a Monte Carlo 
experiment, with time steps of 0.001 and run to a time t = 10 8 . 
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spring constant A, as shown in figure [H So, if we replace / by / = u cos(kx — ut+tp) —u , 
there is a range of uq for which, for sufficiently small or sufficiently large A, the dumbbell 
drifts with the wave, but for intermediate A, the dumbbell drifts against the wave. This 
begs the question of whether there is some mathematically similar flexibility effect in 
other applications, for example in molecular motors. 

For the same reasons as discussed in [2J, these one- dimensional results can be inter- 
preted as components of higher- dimensional results, as at order e 2 , cross-terms for waves 
of either different spatial or temporal frequencies average to zero. Thus the 'fanout' of 
particles of different diffusivities discussed in j2j under the action of several waves would 
extend to dumbbells of the same diffusivities but different spring constants. 

It remains to be seen whether this effect would be significant enough to have practical 
applications. Note, however, the Brownian forcing does not necessarily have to be real 
thermal agitation, but could be some other random forcing, for example eddy diffusion, 
in which case the dumbbell would represent a much bigger flexible body. 
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